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We compare the scattering amplitude resulting from the several quasipotential equations for 
scalar particles. We consider the Blankenbecler-Sugar, Spectator, Thompson, Erkelenz-Holinde and 
Equal-Time equations, which were solved numerically without decomposition into partial waves. 
We analyze both negative-energy state components of the propagators and retardation effects. We 
found that the scattering solutions of the Spectator and the Equal-Time equations are very close 
^^ , to the nonrelativistic solution even at high energies. The overall relativistic effect increases with 

>^ ■ the energy. The width of the band for the relative uncertainty in the real part of the scattering T 

matrix, due to different dynamical equations, is largest for backward-scattering angles where it can 
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O , be as large as 40%. 
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I. INTRODUCTION 



Although hadrons are small bags of confined quarks, at low and intermediate energies what is observed are mesons 
and baryons, suggesting that the genuine QCD degrees of freedom may not be adequate for the study of nuclear 
structure at this energy regime. At low energies the nuclear dynamics is described in terms of nucleons interacting 
via potentials, within the framework of the Schrodinger or Lippmann-Schwinger equations. At intermediate energies, 

f^ , however, the Non Relativistic (NR) quantum mechanical approach may fail. In particular, for processes with high 

[~^ momentum and energy transfer, we may expect the NR scattering equation to be inadequate to describe the nucleon- 

^^ nucleon (NN) scattering. 

^"^ • The NR scattering equation, or Lippmann-Schwinger (LS) equation, gives the NN scattering amplitude T, as: 

-B : r(p', p; W) - F(p', v;W)- jTT^Vip', k; W)- ^— T(k, p; W), (1) 

r- { ' where m is the nucleon mass, p, p' and k are respectively the initial, final and intermediate relative S-momenta; W is 

(— I ' the total energy which depends on p^, and V is the interaction kernel. This kernel, in a NR framework, describes the 

r* 1 instantaneous interaction or potential. The homogeneous term involves also an effective 2-particle scalar propagator. 

. 5^ ] For a relativistic formulation, we have several alternative equations since there is not a unique description dictated 

\^ . from first principles 0]. An implementation of relativity which satisfies covariance explicitly leads to the Bethe- 

$H ' Salpeter (BS) scattering equation M: 



T{p\p; P) = V{p\p; P) + i j ^iV{p', k; P)G{k; P)T{k,p- P). (2) 

Here p, p' and k are respectively the final, initial and intermediate relative 4-monienta, and P the total system 4- 
momentum. As for V , the interaction kernel, it consists of the sum of all irreducible interaction processes derived 
from a considered Lagrangian. In the scalar case the 2-particle propagator G is given by 

^^^' ^^ " [m2 - (P/2 + ky - ie] [m2 - (P/2 - kf - is] ' ^^^ 

Contrarily to the LS equation, which is 3-dimensional and includes only instantaneous interactions, the BS is a 
4-dimensional integral equation and its kernel includes retardation and a dependence on the virtual intermediate 
states energy (off-mass-shell effects). Since the exact solution of the BS equation with the complete kernel is still 
out of reach of present calculation capabilities, several approximations have been developed. Nevertheless, all these 
different methods have in common the truncation of the kernel, which keeps only a restrict set of meson exchange 
diagrams. In Ref. 0, the 4-dimensional BS equations was applied by Tjon and collaborators to the NN problem 
with a kernel derived from One-Boson-Exchange (OBE) generating the so called Ladder Approximation. It has been 



proved, however, that the Ladder Approximation does not verify the one-body hmit Q. The addition to the OBE 
kernel of a crossed-box diagram, implementing a particular off-mass shell extrapolation, has been done by Theufil 
and Desplanques g for a scalar bound state problem. Alternatively, for the same problem Phillips and Wallace Q 
simulated the cross-box diagram through a specific modification of the 2-particle propagator and the result is known 
as the so called 4D Equal-Time equation. Calculations considering ladder and crossed-ladder diagrams have been 
performed by Nieuwcnhuis and Tjon [[7|, based on a Feynman-Schwinger formalism and Monte-Carlo techniques. The 
calculations are, however, still restricted to scalar particles and the bound state cases. 

Even with a truncated kernel the BS equation is, beyond Ladder Approximation, very hard to solve. This difficulty 
has triggered the development of covariant 3-dimensional reduction equations. These can be classified into two classes: 
the QuasiPotential (QP) type equations and the Klein type equations. 

In the QP equations the 3-diincnsional reduction is achieved by fixing the relative energy fcg in ^ covariant way. 
These equations can formally be written as a BS equation, if we replace G by another propagator g including a 
covariant energy constraint. In the last years two QP equations have been applied extensively to hadronic physics 
problems: the Blankenbecler and Sugar (BbS) |q-|lO|] and the Spectator (Sp), or Gross equations |11 Iz 



The BbS equation |13|, is constructed assuming that, in all intermediate states, the two particles are equally off- 
mass-shell. The Sp equation W, is obtained restricting one particle to its positive-energy on-mass-shell state in all the 
intermediate states. Although this equation is exact for scalar particles in the one-body limit, when applied to the NN 
interaction it implies approximations. Furthermore, to account for the identity of the two particles, a symmetrized 
propagator has to be used |11| . In the BbS equation the interaction is instantaneous because the exchange boson do 
not carry any energy, since both nucleons are equally off-mass-shell, whereas in the Sp equation energy exchange is 
allowed and therefore retardation is included. Another QP equation, the 3D Equal-Time (ET) equation, was proposed 
by Mandelzweig and Wallace to deal with the QED problem of spin 1/2 particles [Oj, and with spin zero and spin 
1/2 particles with boson exchange cases |L5| . In this equation the interaction is also instantaneous and includes 
the crossed-box diagram in the eikonal approximation, or high energy limit 0,04-06], where the forward scattering 
components are dominant. This 3D ET equation has been applied by by Tjon and collaborators to the NN system 
with an One-Boson- Exchange kernel |1^ . Calculations have shown that the 3D ET equation results are the closest to 
the BS with OBE and crossed-box kernel for the bound state of scalar particles W. Henceforth we refer the 3D ET 
equation as ET equation. 

The Klein equation, obtained also from the BS equation, is related with time ordered perturbation theory without 
anti-particles g,|l^, and was used by the Nijmegen [|l^,|9| and Bonn ||^ groups in their applications to two- meson- 
exchange calculations. 

In this work we study relativistic effects included in the most well known QP equations - BbS, Sp, Thompson (Th) 
pO| and Erkelenz-Holinde (EH) pl| - and compare their results with ones obtained with the NR Lippman-Schwingcr 
equation. Since the exact solution of the relativistic equation is not known, there is no exact relativistic answer to be 
taken as a reference. Therefore, we focus on the analysis of the relativistic content of the different QP equations. Our 
natural reference, consequently, is the nonrelativistic LS equation. We emphasize, at this point, that the BS equation 
in ladder approximation, which has been taken as a reference in the past |22| , has been shown to be inadequate for 
that purpose because, on the one hand, does not reproduce the one-body limit, as mentioned before, and, on the 
other hand, excludes non negligible contributions of the crossed-box diagrams, as found by the authors of Ref. 0]. 

As a first step, we consider in all our calculations (corresponding to different equations), only scalar particles 
interacting through the Malfliet-Tjon potential [p3| . Only by using the same interaction in all the dynamical cases we 
are able to withdraw conclusions about the way the different formalisms introduce relativistic effects. The relativistic 
effects considered come from kinematics, retardation and negative-energy state components of the propagators. We 
point out that retardation is automatically excluded from the instantaneous BbS equation. As for the negative-energy 
state effects, the different equations account for them in different ways, as it will be shown. 



HI 



In Sec. Q we present the QP formalism, and give the explicit expression of all the propagators. In Sec. 
the method of Ref. |g^ , reviewed in references [|4| , to distinguish between the equations that include retardation from 
the others where the interaction is instantaneous. We also introduce the representation of Wallace and Madelzweig 
p,n5yiqj29] to separate the positive and neg ative-energy state effects, and discuss the suppression of the negative- 
energy states contributions. In reference |22f| a similar derivation was presented, but none of them discuss explicitly 
retardation effects, which is one of the main goals of the present work. In Sec. |^ we explain the numerical method 
used to evaluate the T scattering amplitudes, which does not involve a partial wave decomposition. This method is 
based on the work of Elster and collaborators pq| in a NR framework. In Sec. M we show and discuss the results for 
the scattering amplitude at two different energies. As an example we also show the phase shifts for the s-wave. In 



Sec. VI we summarize and draw some final conclusions and remarks. 



II. QUASIPOTENTIAL EQUATIONS 

The BS equation (^, can be written in a shorthand notation as 

T = V + VGT. (4) 

A common approximation consists on restricting V to the One Boson Exchange (OBE) diagrams, and is usuaUy 
cahed the Ladder Approximation. There are two main objections to this approximation: firstly, as mentioned before, 
the one body hmit is not recovered when one of the masses goes to infinity B; secondly, it has been shown that, in 
general, crossed-ladder diagrams may have important contributions [^ and hence should not be neglected. 

Alternative ways consists on the replacement of BS Eq. (H) by the set of equations: 

T ^ K + KgT (5) 

K = V + V{G - g)K, (6) 

which is perfectly equivalent to the BS equation, and where g is a different 2-particle propagator, containing the 
3-dimensional reduction. Since there are no first principle rules to dictate how this reduction should be performed, 
a large ambiguity gives room to many options. Usually, different approximations compromise different appropriate 
choices of the constraint on the energy-component of the 4-momentum in the propagator g, and different truncations 
of the kernel K. 

The rate of convergence of the new kernel K is determined by the difference G — g. If g is a good approximation 
to G, then only the first term in Eq. (|g) can be kept and Eqs. (||) and (^ reduces to one simpler 3-dimensional, but 
still covariant, equation. This is the standard procedure in any application of QP equations. If needed and possible, 
the lowest order truncation of the kernel K — V can be corrected by higher order terms. 

Any finite order truncation of the kernel beyond the OBE approximation generates unphysical singularities in the 
scattering equation, due to the delta function constraint on the energy J^ . For this reason 3-dimensional QP equations 
with kernels not restricted to the OBE form became controversial P,16l23]. However, we point out that it has been 
shown that the contribution of the two-pion subtracted p7| and crosscd-box diagrams can be simulated by effective 



to and p vector meson exchanges 27 1. This result supports the usefulness of the OBE approximation. 

In the present work we apply several QP equations, with an OBE kernel, to the scattering of scalar particles. More 
specifically, we considered the scalar Malfliet-Tjon potential and took the form of this potential as the kernel of all 
the equations considered. 

We perform the calculations in the cm. reference frame, where the expressions are simpler. In this frame the 
Bethe-Salpeter to 2-particle propagator is written as 

G(k- P) = ^^^ (7) 

^ ' ^ [El- {W/2 + fco)2 - ze] [El - {W/2 ~ k^Y - le]' ^ ' 



Here k represents the magnitude of k, E]^ = \Jm?- -I- k^ and W = \jP^ = 2Ep, if both initial particles are on-mass- 
shell. This expression differs from the original of Eq. (||) by the factor (2?7i)^, which is introduced for convenience. 
In fact, in the scalar case both the kernel and the amplitude T are dimensionless, whereas the amplitude generated 
by the LS equation, which refers to fermions, has dimensions of inverse of energy squared. In order to compare the 
solutions of the two equations, we must redefine V and T of the scalar equation by dividing them by (2?7i)^, meaning 
that we have to redefine the propagator by multiplying it by the same factor. This procedure implies the use of 
dimensionless couplings for both the scalar and Dirac cases. 

As mentioned in the introduction, in the present work we consider the following QP equations: the Sp, Th, BbS, 
EH and the ET. All these equations can be obtained replacing the propagator G of the Eq. (0) by the following g 
functions: 

9Mk;W)=^2.^^^i^ (8) 

E]^W E]^- ^ -le 

,...(fc;I^)=.2.f^45^ (9) 

m md{ko - -Bk + W/2) 



gEH{k; W) ^ z27r- '-^^-^^ ^^ (10) 

\ 4 

m m S{ko - E]^ + W/2) 
KW Ei,-f-ze 



gsp{k;W)=i2Ti— — — ^ -—. (11) 



The BbS and EH propagators are identical except for the argument of the delta function containing the energy 
constraint. The same can be said about the Th and the Sp propagators. We also consider the ET scalar equation 
|P5[ , corresponding to the propagator 

,Mk;W)=^2.-^^-^^(2-^J, (12) 

that differs from BbS propagator (g) by a kincmatical factor, which becomes one when the particles are on-mass-shell 
{W = 2E^). 

In all cases the energy-fixing condition is included through a (5-function. We note here that even the NR equation 
(|l|) can be obtained from the BS Eq. (0) by replacing the 2-particle G propagator by the NR propagator 

9Nnik■,W)=^2n^,!I^^^^, (13) 

which can also be written as 

g^n{k■,W)=^2n J^f:>^ . , (14) 

but this does not mean that the LS equation is covariant. 

All the propagators defined above are for scalar particles. The first five propagators can be organized into two 
classes, by means of a function f{k;W), the BbS and Sp types ||l|, p^ , p8[ . The first class includes instantaneous 
equations while the second one includes equations with retardation. 

gB,sik;W)^^2n^-^^i^^^6iko) (15) 

gsp{k; W) ^ ^27^^ rT^^w^\.J iko -E^ + W/2). (16) 

Any choice of / is possible providing that 

(17) 

which results from the simultaneous constraint of the 2 nucleons in its physical state (positive-energy on-mass-shell 
state). Taking /(k; W) = 1 in Eqs. ( [l5| ) and ([l6|), we obtain respectively the original BbS and EH propagators. Taking 

mw) = E^, (18) 

in Eqs. (|l5| ) and (16), we get respectively the Th and Sp propagators. Finally the ET propagator is obtained from 
Eq. (p^) using the following f{k]W) function: 



E^El-^~^e 


m mf{k;W) 


£^k El X '^ 


f{p;W) = l, 



w 



/(k;W^)=(^2-^j. (19) 

We point out that the function /(k; W), to a great extent arbitrary, is related with the weight of the contribution 
of the negative-energy state components, as we shall see. 

For numerical applications we rewrite the propagators (||)-(0) and (n2h in a different form. Firstly, we represent 
by g the propagator without the (5-function factor. Secondly, we define the nonsingular function /(k; W) accordingly 
to 



f{k;W) 

'" = E 

and /(k; W) is expressed in terms of /(k; W) by 



-g{k;W)= ^^ '^ ' (20) 



mm^fOc^^ (21) 



These /(k; W) functions verify the relativistic elastic cut condition, expressed in Eq. (|l7|), corresponding to 

mW) = 2{^y. (22) 



III. NEGATIVE-ENERGY COMPONENTS OF THE PROPAGATORS AND RETARDATION EFFECTS 

IN QP EQUATIONS 

In order to clarify the interpretation of the positive and negative-energy state effects we follow here the analysis of 
Ref. [Q. We start with the homogeneous term of the Eq. (H), in the cm. frame where P = {W, 0) 

7^(k; P)= j j d^kdfcoVb^, p'; fco, k; W)G{k- P)T{ko, k;po, P; W); (23) 

in the above equation we explicitly separate the energy and 3-momentum components of all momenta. Using the 
representation of Wallace and Mandelzweig |H,Eq,lq|, the propagator of 1-particle with 4-momcntum P/2 + k, can be 
written as 

G, (k- P) = - f ^1^ ^r ] .24. 

'^ '^ \ko-E^ + W/2 + ie ko + E^ + W/2-iej' ^ ' 

and the propagator of particle 2 with of momentum P/2 — fc as 

G,(fc:P) = ( ^ ^ 1 . (25) 

-^ ■ ' \kQ-Ei, + W/2-ie ka-Ei,-W/2 + ie} ^ ^ 

For scalar particles N^ = N^ — ra/ E]^ whereas for fermions N/^ and N~ are respectively the positive and negative- 
energy projectors |6 p^ , p^ . 
Defining 

W 
uj+^E^- — ^ie (26) 

W\ 
E^ + —j+ie. (27) 

we can rewrite the BS 2-particle propagator G as 



'< fcn — OJ, kc\ — bJ, > 



Ei>0 Ei<0 

^ ko+uj^ ko+Ljy_ i 

B2>0 B2<0 

In the expression above we indicate, for the two terms of each particle propagator, the sign of the particle energy at 
the corresponding pole. 

Following [ p2|p4 |, the QP equations can be obtained in a different way than the one discussed in the previous 
section, namely by fixing the relative energy fco in V and, consequently in T . If we choose fco = in U and T, we can 
factorize the energy and B-momentum integrations in the following form 

7^(k; P)^ j d3kU(0, p'; 0, k; W^)r(0, k; 0, p; W)5(k; P), (29) 

where Q is 

gik;P)^ fdkoGik;P) (30) 

With this choice the interaction between the nucleons is instantaneous. We call these equations the instantaneous- type 
equations. The integration in the fco variable can be done analytically using the residues theorem . The pole structure 
of the integrand function is presented in the Fig. nl Choosing the lower half contour for the fco integration, we can 
conclude that 

g{k;P) = h{u;+) + l2{-u;n, (31) 



where the residues Ji and I2 are 

^1^^ \ W W + 2E^} ^ ' 

By keeping both Ii{iL!^) and /2(— wr) in Eq. (£9) we obtain the BbS equation. Neglecting the residue l2{—uj7) we 
get the Th equation. 

The other class of equations can be obtained when we fix the fco in order to correspond to the positive-energy state 

4- 



of the particle 1. This corresponds to the choice fcg = uj^ . So, in that case 



7^(k; P)= J d^kV{u;+,p'; u;+, k; W)T{u;+, k; u;+ , p; W)g{k- P). (34) 

All equations in this class include retardation, because they includes the transfered energy term in the interaction 
kernel. If in Eq. ( p4|) we include the two residues Ii{uj^) and l2{—uj^) we obtain Erkelenz-Holinde (EH) equation, 
and if we neglect the last one, we obtain the Sp equation. In that respect Sp equation is more consistent since it 
considers the same value of fco in all integrand functions. 

We summarize in Table m the comparison among the several QP propagators accordingly to the the discussion 
above. 

Using the decompositions (^2|) and ( |33| ) we may write 



„ , „ , r N+N+ N+N- NZN+ N~N^ ] , , 

g...(k; W) = fe„(k; W) = ,:2, I ^^^ ^ ,;_ _^ + ^^ - ^^ + ^^X_l_| (35) 

To avoid any possible confusion we call the reader's attention to the fact that, for spin 1/2 particles, some authors 
call BbS to a somewhat different 2-particle propagator |ll|,|6|,p|jl3[ . 

For the BbS/EH equations we need to consider the sum of both terms in (pTl). It is interesting to note that, if the 
weight functions are equal, N^ = N~ , as it happens in the scalar case, the residue /i at the negative-energy pole of 
particle 2 cancels exactly against the residue I2 at the positive-energy pole of the particle 1. As a result, we can write 
for the scalar case 

gBbsi^lW) ^ i2TT ( —] I ^ + ^ 1. (37) 

^ ^^ ' ' \E],J \2E],-W-ie 2E], + W j ^ ' 

where -q — 1 corresponds to the maximal inclusion of negative-energy state term of the propagators, and the case 
77 = totally excludes them. In the same situation we obtain for the Sp/Th equations the following 

fc,(k;«')^.2,(£)'|5^^-l^ + ^}. (38) 

Comparing (|37| ) with (|3^) in the case of negative-energy suppression (77 — 0) we conclude that all the propagators are 
equal. 

As for the ET propagator, it has been obtained in Refs. Mjlq] and reads 



ferik; W) = .2. 1 ^^^ ^^f_^^ ^^-^^^-^^ _^ I . (39) 

Compared to the BbS propagator, the qet contains two additional contributions, which the authors justify as an 
effective inclusion of the crossed-box diagram. Restricting to the scalar case, this expression simplifies to 

,,H.<;H')^.2.(|i)^{,5^3Wu + i+2l^}^ (40) 

If we restrict ourselves to the positive-energy components (77 = 0) then the ET propagator is equivalent to all the other 
propagators. In particular, for a complete suppression of negative-energy states, BbS and ET equations generate the 
same results. 
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FIG. 1. Poles structure of the G propagator. 
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TABLE I. Comparison among QP equations. Here fco 
suppression of the residue. 



LJ^ means that retardation is taken into account, and x means 



IV. SOLVING THE SCATTERING EQUATION WITHOUT PARTIAL WAVE DECOMPOSITION 

In the previous section we derived the QP equations in a way to make exphcit, in each of the considered equations, 
the contribution of the positive and negative-energy terms of the BS propagator, as well as the effects of retardation 
present in the kernel. However, as mentioned before, this method generates exactly the same QP equations as the 
one described in Sec. O, and we chose to proceed here with the formalism presented in the latter section. 

Our starting point is, consequently, Eq. (j^) with g given by one of the functions of Eqs. (p|)-([l2|), which after 
performing the /cq integration can be written as: 



T(p', p; W) = F(p', p;W)^ J A^y(p', k; W)g{k; W)T{k, p; W). 



(41) 



Here g is given by Eq. (|2C|). The value of the relative energy variable, which is not explicit in Eq. (|4l|), is fco = for 
the instantaneous-type equations, and kg — E]^ — W/2 for the retarded-type equations. 

The QP equations are 3-dimensional integral equations, corresponding to a radial integration and a 2-dimensional 
angular integration. Since the physics is independent of the scattering plane, we may choose the incoming 3-momentum 
p to have the z-direction and the scattering plane as the xz plane, meaning that the final momentum is only defined 
by its magnitude and by the angle 9, which then gives the scattering angle. If the interaction is scalar (spinless 
particles), the T- matrix and the kernel V are independent of the angle ip that specifies the scattering plane. 

Representing the magnitudes of the initial and final 3- momenta by p and p' respectively, and defining u = cosO, 
we can write then 



T{p',p;W)=T{p',u;p,l;W), 
V{p',p;W)^V{p\u;p,l;W). 



(42) 
(43) 



The intermediate 3-momentum is characterized in terms of a polar angle, v = cos^j., and an azimuthal angle ip]^. 
With this notation we write 



T{k,p;W)^T{k,v,p,l;W), 
V{p\k;W)=V{p\u;k,v,^l,;W). 

Using Eqs. ( [43|) and (^) and factorizing the integration in (pi^ we can rewrite { \U\ ) as 

T{p',u;p,l;W) = Vip',u;p,l;W) 
k^dk '•1 



(2n) 



dv V{p', u; k, v; W) ■'^^^'^^ . r(k, v; p, 1; W) 



(44) 
(45) 



(46) 



where 



V{p',u;k,v;W) 



1 



V{p',u;k,v,(p]^;W)d(p]^. 



Treating the singularity of the integral function in the usual way, we write 

T{p\u;p,l;W) ^V{p',u;p,l;W) 



, _ f^ k^dk 

dvV -—. 

-1 Jo (Stt, 

dv 



(2rn)2p 
W 



i(27r) 



^V{p',u;k,v;W){^^^T{k,v;p,l;W) 
Vip',u;p,v;W)T{p,v;p,l;W) 



(47) 



(48) 



where we used 



kV"(k;P) 



d /p W] 
TTl? l^k - -J 



(2m)2p 



k=p 



(49) 



accordingly to the eq. (p2|). 



By standard discretization we transform the integral equation into an algebraic linear set of equations, as explained 
in the appendix |A| and obtain, for a fixed on-mass-shell initial state, the half-off-niass-shell scattering amplitude 
T(p^,u,;p,l;W^). 

One of the numerical checks that has been performed is the verification of the optical theorem. The number of grid 
points for the magnitudes of the 3-niomenta and for the cosine of the polar angles are fixed by imposing an accuracy 
of the order of 1% in the verification of the optical theorem. Typical values used were 17-25 for the magnitude of 
the momenta and 20-50 for the angular discretization. We point out however, that the convergence of the on-shell T 
matrix does not require so many grid points. 

V. RESULTS AND DISCUSSION 

We solve the scattering equation (||) with the propagators (|q)-(|l^). For instantaneous- type equations we considered 
the kernel to be the Malflict-Tjon potential p3]: 



Vr V^ 



^(k, P; W) ^ „. " ., - ,. , ,., „,, , (50) 



with the parameters of Ref . 1 26 



yR = 3,1769 /ij? = 305, 86 Mel/ (51) 

Va = 7,9210 A*A = 613,69 Mey (52) 

The choice of the interaction is due, on one hand, to its simplicity, and, on the other, to the fact that it as has been 
used recently in the literature pq] , within a nonrelativistic framework and without a partial wave decomposition for 
the T-matrix. This way we could compare and control our numerical results. 

In order to study retardation effects however, we have to consider a covariant version of Eq. (pQ), obtained by 
replacing the square of the 3-momentum by the negative square of the 4-momentum (q^ -^ ^'f')- Therefore we use 
for retarded-type QP equations 



Vr Va_ 

1^1 + {\^- p)2 - {E^ - E^Y A^i + (k - p)2 - [E^ - E^y 



V{Kp;W)= .., , „, „,,^ ,^ .TTI - .,2 , n. „^. .z. ETTI ■ (53) 



To evaluate the scattering amplitude without partial wave decomposition we apply the method presented in the 
appendix^. We consider throughout our calculations two energies: Tiab — 300 MeV, where we expect small relativistic 
effects, and Tiab = 1 GeV, which corresponds to the intermediate energy range. The results are presented in the Figs. 
g, g and |j. We will focus mainly on the real part of T, where the effects are more significant, but the imaginary part 
is also presented for completeness. 

We organize our study accordingly to the following topics: 

• Effects of the negative-energy components of the propagators- where we compare the QP amplitudes 
calculated with positive-energy components of the intermediate propagators [i] ~ 0), with the corresponding 
amplitudes including the full propagators {rj = 1). Moreover, by comparing the amplitudes EH/BbS with the 
Sp/Th ones, we separate the effects originated from I2{—ujZ), present only in the first pair of equations. We 
note that accordingly to Eqs. ( p^ ) and (|3^) the inclusion of the I2 residue means extra terms originated by 
negative-energy components. In particular, only I2 contains a term involving negative-energy components of 
the propagators of both particles. 

• Effects of retardation- where we compare the results of Sp/EH equations with and without retardation. We 
note that, as shown above, the Th equation is the instantaneous version of the Sp equation, and that BbS 
equation the instantaneous version of the EH equation (Table I). 

In Fig. 0we show the Sp and Th scattering amplitudes calculated with (77 = 1) and without (77 = 0) negative-energy 
components of the propagators. The + label corresponds to 77 = 0. As can be seen from the figure, the Sp and NR 
amplitudes are indistinguishable at 300 MeV and very close at 1 GeV. Since the Sp equation includes the combined 
effects of retardation and negative-energy state components, we may conclude that the two types of effects tend to add 
and to approximate the Sp relativistic calculation to the NR one. We can also see that the negative-energy components 
decrease the real part of the scattering amplitude for Sp and Th formalisms, and, in the former, retardation enhances 
this effect. These conclusions can be drawn for the two energies considered. We point out that the difference between 
the Sp and NR amplitudes increase with energy as expected. 



In Fig. ^ we compare the EH and BbS scattering amplitudes, which differ only through the inclusion of retardation 
in the former. It is also shown the ET result. In all cases the inclusion of negative-energy components decreases 
the real part of the scattering amplitude, and approaches the amplitudes to the NR limit. More significative is 
the ET result, which has more content of negative-energy components than the EH one, and is closer to the NR 
result. We remind here that the ET propagator without negative-energy states corresponds to BbS+ (BbS with only 
positive-energy states). By comparing the retardation effects that are present in EH and absent in BbS, we conclude 
that retardation also decreases the real part of the amplitude. This is valid for the positive-energy versions and for 
the full versions. Once again, we may conclude that the combined effects of retardation and negative-energy state 
components tend to approximate the relativistic calculations to the NR ones. The fact that the ET curve, which does 
not contain retardation, is the closest to the NR result is explained by the larger contribution of the negative-energy 
state components in its propagator. 

The full scattering amplitude results for the five QP equations along with the NR one can be observed in Fig. ^. The 
curves shown confirm the conclusion drawn above, namely that the combined effect of retardation and negative-energy 
state contributions tend to approximate the relativistic calculations to the NR ones. We point out, however, that the 
full result depends upon a delicate interplay of these two relativistic effects, and, in particular, on the amount of the 
contribution of negative-energy state components of the propagators. This can clearly be seen by comparing the Sp 
result (no l2{—uj^)) with the EH result (with l2{—uj^)), being the first the closest to the NR limit. 

All the previous considerations about the real part of the scattering amplitude remain valid for the imaginary part 
with an opposite direction that is, in general retardation and negative-energy effects increases the imaginary part of 
the scattering amplitude. 

We note that the result from the BbS equation which, accordingly to Eqs. (^) and (|lj), indeed generates the NR 
equation when we replace -g^ by 1, deviates, the most from the NR result. 
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FIG. 2. Spectator (Sp) and Thompson (Th) scattering amplitudes for 300 MeV and 1 GeV. NR corresponds to 
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FIG. 4. Comparation between all amplitudes, for 300 MeV and 1 GeV. 

Finally, as an example, we show in Fig. H the phase shifts for the I — partial wave, generated by the different QP 
equations, as a function of the incoming energy. As can be seen from the figure, at low energies the different equations 
generate similar results, but at higher energies discrepancies are obvious. In this partial wave the phases generated 
by the EH and Th equations are the closest to the NR ones, contrarily to what happens with the full T- matrix, where 
the Sp and the ET are the closest. This fact has a very simple interpretation, stemming from the fact that at high 
energies the contribution to the transition matrix of the more peripherical partial waves becomes quite important, 
and therefore looking only to the s-wave channel is a very limited analysis. This conclusion is also achieved in the 
workof Ref. M. 
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FIG. 5. s-wave phase shifts generated with the different QP equations, as a function of the energy. 

VI. CONCLUSIONS 

We consider a family of relativistic quasipotential scattering equations, and study the effects of retardation and 
negative-energy state contributions in the propagators. We restrict our calculation, for the time being, to the scalar 
case. Comparison with the NR equation is provided. The numerical calculations of the T matrices were done without 
partial wave decomposition. Instead, we solved the two-dimensional integral equation, and the accuracy of the method 
was checked through identity provided by the optical theorem. The method implied a discretization of the integral 
variables with 17-25 grid points for the magnitude of the relative momentum, and 20-50 grid points for the angular 
variables. The NR solutions were compared with the results of Ref . |g^ , where also the partial wave decomposition 
has been avoided. This comparison provided an extra check to our method. 

We conclude that, in general, retardation decreases the real part of the scattering amplitude, and increases the 
imaginary part of the same amplitude. 

The inclusion of negative-energy state contributions induces a similar effect. 

The equations that include the I2 term, i.e. the ones containing additional negative-energy state effects, deviate 
the most from the NR results. We found also that the Sp and the ET equations generate the scattering amplitudes 
which are the closest to the NR result. We also point out that the Sp (with retardation but no I2 terms) and the 
ET equations (without retardation but with I2 terms) consider, in an effective way, some contribution from the 
crossed-box diagrams. On the other hand, the ET equation, on top of the I2 term, includes further negative-energy 
state effects, associated with its specific way of describing effectively the crossed-box diagrams. These additional 
terms compensate partially the absence of retardation. We may therefore conclude that the inclusion of crossed-box 
diagrams, even in an effective way, approaches the QP results to NR ones. 

It is interesting to note that the way the Spectator formalism deals with negative-energy state components, retar- 
dation and effective crossed-box diagrams, makes the corresponding amplitudes the closest to the NR limit, even at 
high energies. Finally we note that the relative uncertainty in the real part of the scattering T matrix due to different 
dynamical equations depends on the scattering angle and is the narrowest near the forward-scattering region. In the 
backward-scattering, it can be as large as 40%. This could be expected since in this region the momentum transfer is 
larger. 

In this work we restrict ourselves to spin zero particles. However, calculations to account relativistically for the 
spin 1/2, such that our formalism can be applied to the nucleon-nucleon system, is on progress. In such a more 
realistic calculation, the uncertainty band found in this work for the T matrix will most likely imply differences 
in the nucleon-nucleon interaction, compatible with the scattering data and originated by the different dynamical 
formalisms. 
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APPENDIX A: SOLVING SCATTERING EQUATION WITHOUT PARTIAL WAVE DECOMPOSITION 



In order to solve the 2-dimensional integral Eq. (^8|) in the variables k and v, we discretize this variables using a 
gaussian grid with the following properties: 

• Linear grid of N^ points Uj the variable u G ] — 1, 1[, with weight hj. 

• Linear grid of Np even points Xi the variable x £ ]0, 1[, with weight Wi. The correspondent momentum points 
are obtained through the transformation 

p,=A^, (Al) 

"■•' - ^ ^ w., (A2) 



where A is a free parameter. We and add the point XNp+i — 1/2 with w'j^ +i = 



We take A — p the on- mass-shell initial momentum. With this choice the evaluation of the principal part in (48) is 
simplified 

N„ 

(A3) 



pOQ P 

V /(k)dkc.^/(p,K 

•^0 i^i 



The point pjy +i = p corresponding to the pole is not included in the sum, hence, the function have no singularities. 
We take the outcoming variables p' and u from the same grid. In that case the integral equation (Esh is transformed 
into a linear algebraic system of equations with the unknown ^(y) = T(pj,Uj;p, 1;M^). The brackets in ij indicates 
that this indices can be contracted in only one. The linear system can be written as 



where 



The matrix C are given by 



where / in the unity and 



C ■T = V, (A4) 

%)=y(p„Wj;P,l;M^). (A5) 

C = I + A + iB, (A6) 



{ITT) hp^, 2" 

1 (2to)^p - 

^(v)Xi'r) = ^ ^ V{P^,Uj;P^',Ujr,W)6i,.Np + l (A8) 

The factor Si'^N +i gives the elastic-cut constraint and fixes the momentum Pj, = p = p^ _^i. Note that A is zero if 
i' — Np + 1, and B is always zero except for i' — Np + 1. Accordingly, the matrix elements oi A + iB are purely real 
or imaginary. 
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